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Abstract 

A mathematical model describing motion of an inhomogeneous incompressible 
fluid in a Hele-Shaw cell is considered. Linear stability analysis of shear flow class 
is provided. The role of inertia, linear friction and impermeable boundaries in 
Kelvin-Helmholtz instability development is studied. Hierarchy of simplified one¬ 
dimensional models of viscosity- and density-stratified flows is obtained in long-wave 
approximation. Interpretation of Saffman-Taylor instability development is given 
in the framework of these models. 
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1 Introduction 

A general property of fluid shear flows is the Kelvin-Helmholtz instability. This type 
of the interface instability of ideal fluid layers moving with different velocities is studied 
clearly enough. Classical examples of linear analysis of such flows are described in DEI 
and in others monographs. Recent review [3] is devoted to stability analysis of viscosity- 
stratihed shear flows. A substantial theoretical object having practical applications is 
the Hele-Shaw shear flows [1]. Displacing more viscous fluid by less viscous one leads 
to formation of viscous hngers. It is caused by Suffman-Taylor instability [5]. Classical 
models of fluid motions in a porous medium are based on linear Darcy law However 

shear instability can develop on the longitudinal interface in the case of fast flow at formed 
viscous hngers also [5]. It is necessary to use the nonlinear Darcy law for analysis of this 
phenomenon. It takes into account inertia terms in the momentum equations. Such 
models are considered in [91 HD]. There are exact and experimental studies of viscous 
hngers formation with inertia ehect for immiscible huids in the papers cited. The same 
problem for miscible huids is described in m- It is shown that inertia slows down the 
growth of viscous hngers. Ehect of Hele-Shaw cell constriction and walls elasticity on 
Sahman-Taylor instability is analyzed in nans]. 
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In the manuscript the nonlinear system of equations describing inhomogeneous fluid 
flow in a Hele-Shaw cell is under study. The model admits a class of shear flows. The 
class is characterized by arbitrary dependence of the horizontal velocity component and 
density on the vertical spatial coordinate. One of the main purpose of the study is stability 
analysis of this class of flow. The equation for amplitudes of small perturbations is 
obtained. The conditions on the interface are formulated for fluids with different physical 
properties. Stability analysis of two-layer Hele-Shaw flow is carried out in details. It is 
shown that the interface of layers moving with different velocities is instable for short-wave 
perturbations and stable for long-wave ones. 

The pressure changes weakly in the vertical direction for prevailing horizontal fluid 
motion ng. It allows to consider long-wave approximation and class of layered flows. 
A hierarchy of one-dimension models is constructed basing on two simplihcations of the 
momentum equations. They are linearization and using of Darcy law. Two-layer flow 
of viscosity- and density-stratihed fluid is studied. Additionally the evolutionary form 
of the governing equations is obtained both in multi-layer and two-layer cases. The 
hyperbolicity property is proved for the models above. It enables to apply Godunov 
method modihcations for numerical calculations. In particular Nessyahu-Tadmor central 
scheme is used in the paper. Instability of the interface is proved at displacing more 
viscous fluid by less viscous one if density stratihcation is neglected. 

The paper is organized as follows. In Section 2 we introduce the mathematical model 
of the process described above. In Section 3 stability analysis of shear flows is carried 
out. The cases of homogeneous and stratihed flows are studied. In Section 4 we consider 
layered flows in long-wave approximation. We describe two simplihed models of such flows 
which are treated both analytically and numerically. 

2 Mathematical model 

Equations of motion of a viscous incompressible fluid in a Hele-Shaw cell (region between 
two parallel plates separated by a small gap) have the form 

p(vt-h (v • V)v)-F Vp =-F pg, V-v = 0. (1) 

Here x = {x,y,z) are the coordinate vector, t is the time, v = {u,v,w) is the velocity 
vector, p is the pressure, p is the density, p is the viscosity and g = (0, —p, 0) is the 
acceleration of gravity vector. The operator V is calculated with respect to the coordinate 
vector. The characteristic sizes of cell {L,H) in the x and y direction respectively are 
signihcantly higher than the cell gap b. That is why the summands and \yy vanish in 
the momentum equations. They are negligible compared to v^j,. We assume the density 
and viscosity are monotonic dependencies of concentrations x: p = p(x) and p = p{>c). 
The variable x varies from zero to one and satishes the equation 

Xf -|- V ■ Vx = 0. (2) 
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It should be noted that x can change continuously or take piecewise constant values. 
The hrst case corresponds to an inhomogeneous fluid, the second one characterizes a 
multicomponent fluid. 

We consider the velocity held in the form 

u = ^(l- ^u{t, x,y), n = ^ (^1 - j x,y), w = 0. 

It provides the fulhllment of no-slip conditions on the cell walls 2 ; = iLb/2. We also 
suppose that the functions p and x do not depend on z. Integrating equations ([I]) and 
([2]) from —6/2 to 6/2 leads to the system 

p{ut + (3 {uu^ + vuy)) +Px = -pu, 

p{vt + (3{uvx + vvy)) +Py = -pv - pg, (3) 

Ux + Vy = 0, X^ -|- MXj, -|- V>Cy = 0. 

The primes are omitted. Here and below p denotes the modihed huid viscosity 12p/6^. 
The coefficient f3 is equal to 6/5. The simplihcations above are connected with inequali¬ 
ties b << L, b << H and often used for modeling of Hele-Shaw hows. For the constant 
density and using of linear Darcy law equations ([2]) reduce to the classical problem with¬ 
out dihusion consideration [7]. In the framework of the model described the interface 
instability is provided if displaced huid is more viscous than the displacing one. 

The following Sections are devoted to stability analysis of steady shear solutions of 
equations ([3]). Also a class of layered hows is considered in long-wave approximation. 


3 Stability analysis of shear flows 

Equations ([3]) admit the class of solutions 

u = U{y), u = 0, p = P{y) - ax, x = C{y), 


(4) 


where a > 0 is constant, U, P and C are arbitrary smooth functions. The following 
relations are valid 

d 

^ = = p = p{C)=R{y), P'{y) =-gR{y). 

The vector of velocity components, pressure and concentration of perturbed how deviate 
from the basis how weakly. That is why these values take the form 

u = U{ii)-\-u, V = V, p = P{y) — ax + p, K = C{y) + k. (5) 

Here u, v, p and x are the functions of all independent variables. Dependencies p = /i(x) 
and p = p(x) have the form 

U'{y)k X ^ , R'{y)k 


a 


p = 


(7(S)1 U(y)C'(y))' '' + C 


(») ■ 
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Substitution of solution ([5]) into equations ([3]) and linearization of the system obtained 
give the equations for small perturbations hnding 


n/TT rr/ \ 1 “ fU'x \ 

ut + (3{Uu:, + U v) + -p^ = 

Ux + Vy = Xt + UXx + vC = 0. 


( 6 ) 


The tildes are omitted. The prime denotes the derivative with respect to y. 

We use the stream function x, y) for the perturbations of the basic flow. It satishes 
the equalities u = 'ipy and v = —'tpx- We differentiate the hrst equation in ([2D with respect 
to y and the second one with respect to x. Then we eliminate the pressure p and obtain 
the system for perturbations of the stream function and the concentration only 


(9, + /3(/0.)(a^ + ^4,,) - ((/" + 


aU' /U'x 



U'C" 

C' 


X 

O 


U' 


{dt + f/9x)x - C'Vx = 0. 


(7) 


A solution of these equations is sought in the form 


(V', x) = (T(|/), M{y)) exp {ik{x - cf)), 


( 8 ) 


where c is the phase velocity, A; > 0 is the wave number, 27r/A; is the wave length and i is the 
imaginary unit. Substitution of formulas (|HD into (ED leads to equations for perturbation 
amplitudes T(?/) and M{y). Due to the second equation in ([7D we have M = (D—c)“^C"T. 
The hrst equation in (ED reduces to the second order ODE for the function \1/(|/) 


(^/3U -c 
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kUR 
iaU' 2U — c 


kRU'^ U-c 
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(9) 


If the how region is restricted by the horizontal walls y = ±h then the boundary conditions 
for equation (|9D have the form \l/(—h) = \l/(h) = 0 and correspond to impermeability 
conditions on the rigid walls. In the absence of walls the limit relations T —)■ 0 fulhll at 
y —±oo. 


3.1 Conditions on the interface 

We consider shear how (jlD with discontinuity at y = 0. We do not take into account the 
surface tension. The limit values and R^ of functions U{y) and R{y) are diherent at 
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y —)■ ±0. The superscript “+” corresponds to the upper fluid. The linearized kinematic 
condition takes the form 

[dt + U^d^)C = v^ = 

where C{t,x) is the small perturbation of the interface. Eliminating the function C{t,x) 
from the latter equation gives the formula {dt + U^dx)ilJ^ = {dt + U~dx)'ip^. For wave 
solutions dH]) this formula is rewritten as 

(t/+-c)T-= (t/--c)T+. (10) 


On the interface the linearized dynamic condition has the form 


[p - gRQt = 0, 


where [/]! = /+ — /_ is the difference between limit values at the discontinuous. In order 
to obtain the dynamic condition for the stream function we apply the differential operator 
{dt + U^dx)dx to the latter equation. Taking into account equations ([S]) the kinematic 
condition on the interface and form of solution ([8]) we have the formula 




U+ -c 


M!- (11) 


The solution of equation ([2]) and conditions (ITU]) . fITT]) allow to obtain the dispersion 
relation c = c{k). Analysis of this relation gives the conclusion about the flow stability. 
It should be noted that at a = 0, /3 = 1 formulas (jH]), flTOj) and flTT]l reduce to the classical 
relations of linear stability theory of plane-parallel shear flows of ideal fluid mi- 


3.2 Stability analysis of two-layer flow 

We consider two-layer motion of fluid with the interface ?/ = 0 between impermeable walls 
y = ±h. The flow scheme is presented in Fig. [T] The fluid layers move with constant 
positive velocities Ui, U 2 , have constant viscosities pj and densities pj {j = 1,2). The 
compatibility condition UjPj = a follows from formulas (j4]) and is assumed to be valid. In 
the case of piecewise constant distribution of velocity and density equation ([H]) takes the 
form T" — = 0. We solve this equation with conditions \l/(±/i) = 0 and use formulas 

(ffoD, dm). As the result we have the dispersion relation 

({(^RJl - c)kpi - iaf/f A {Ui - c) -1- {{PU^ - c)kp 2 - 

( 12 ) 

-iaU^^) {U 2 - c) - g{p 2 - pi)p{h, k) = 0, 

where ip{h,k) = tanh(2hfc). When the impermeable walls are absent, i.e. ^(y) —)■ 0 
at 1 / —)■ ± 00 , the dispersion relation has also form fim but (p{h,k) = 1. Formula (im 
connects the complex phase velocity c = Cr + iCi and wave number k at given values of f/j , 
Pj, a and (3. If there are such values of parameter k for which equation 0121) has complex 
roots with positive imaginary part Cj then the considered flow is instable. 
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Figure 1: Two-layers flow scheme. 



Figure 2: Dispersion curves (C/i/172 = 6/5, a = 
ll2): 1 — /3 = 1- 2— P = 6/5. 


3.3 Two-layer flow of a homogeneous fluid 

We consider two-layer fluid with identical densities pi = p 2 = 1- Then the gravity force 
and walls do not influence on the stability of the flow studied. Following [1] the coefficient 
P = 6/5 is often considered to be equal to one. When a = 0 {a = piUi = P 2 U 2 ) 
and = 1 we have q = Imc = ±(17i — U2)/2 from dispersion relation flT^ . These 
dependencies are presented by dashed lines in Fig. [2] and correspond to instability of 
the contact discontinuity in homogeneous ideal fluid. When a > 0 the roots of quadratic 
equation (IT^ depend on the wave number k. Special dependencies c, = Cj(fc) are presented 
in Fig. [2] (lines 1). The plots are obtained at parameters Ui = 6/5, U 2 = 1, a = 1/2 and 
Pj = a/Uj. The values are pointed out for the CGS system. Also they can be assumed 
to be dimensionless. It is justihed because we are interested in qualitative behavior of 
solutions. The upper branch of the graph Cj = Ci{k) increases monotonically from zero to 
{Ui — 112 ) 12 . The long waves {k —)■ 0) are stable neutrally while perturbations of short 
waves grow with the same velocity as at a = 0. 

When /3 = 6/5 the dispersion curves c, = Ci{k) have the same qualitative behavior as 
at /3 = 1 (lines 1 in Fig. [2]) if the following inequality fulhlls 

/ = (1 -f p)^iUi + U2)^ - mui + f/|) < 0. 

At /3 = 1 this inequality holds because —{Ui — t/ 2 )^ < 0. One can readily see that at 
> 1 and considerably small difference \Ui — 1 / 2 1 the inverse inequality / > 0 is valid. 
The roots of equation ffT2|) take the form 

Cl,2 = ^((l + /5)(f/l + 172)±/7) 

in the case of a = 0. It corresponds to instability of two-layer ideal flow with correction 
factor > 1 in the momentum equations. Ata>0, /3 = 6/5 and / > 0 the inherent 
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velocity behavior is presented in Fig. [2] (lines 2). In this case the greatest growth of 
disturbances occurs in the middle range of wavelengths. 

It is interesting that at /5 = 0 the flow is stable. Applying Vieta theorem to the 
quadratic equation 


t/l + 172 


za 


la 

~k 


= 0 


kUiU2^ 

proves that there are no complex roots with positive imaginary part. Actually let Ci = 
Cir + icii and C 2 = C 2 r + ic 2 i be the roots of the latter equation. Then from Vieta theorem 
we conclude the fulhllment of the following inequalities 


I^2r ^Ir !^2i ^ O5 ^\r I^li ^2i/^2r ^ 0* 

Due to the inequalities cir + C 2 r = {Ui + V2)/2 > 0, cir > 0 and/or C 2 r > 0 we infer that 
cii < 0 and C 2 i < 0. 


3.4 Two-layer flow of a stratifled fluid 

We consider the shear flow (Fig. [T]) with density stratihcation p 2 = 1, A = pi/p 2 < 1- The 
function ip{h, k) varies in the interval (0,1) and makes adjustments to the last summand of 
the left part in 0121) . It does not influence on the qualitative results. At a = 0 dispersion 
relation flT^ has the following roots 

- (A + l)(Aq + + ^(1 - A^)). 

As in the previous example there is difference between the cases /3 = 1 and (3 = 6/5. At 
(3 = 1 and a = 0 the flow is instable always because at k ^ oo the following relation is 
valid 

Imc ^ T{Ui - f/2)(l + A)-1\/A. 

It is pointed out in Fig. [21a by dashed lines. At /3 = 6/5 there are parameters for stable 
flow. 

Fig. El a presents the dispersion curves obtained from solution of equation flT^ for 
the parameters (3 = 1.2, Ui = 1.5, U 2 = 1, g = 980, h = 10 and A = 0.99. Lines 1 
correspond to a = 0 and lines 2 are pictured for a = 9. The stability area is essentially 
wider compared to the homogeneous flow because of the positive values Cj only for /c > /c*, 
where k^ is the critical value. In the example under study fc* ~ 83 at a = 0 and /c* ~ 55 
for a ^ 0. The values k^ increase monotonically with the growth of the shift value 
\Ui — 1 / 2 1. When the parameter a decreases curves 2 draw close to lines 1. We note that 
at (3 = 1 the qualitative behavior of dispersion curves is the same. The values of critical 
wave numbers coincide for a = 0 and a ^ 0. For the example considered these values are 
fc* ~ 79 and 76. We take Ui = 6/5 and the other parameters do not change. This case 
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Figure 3: Dispersion curves: a — a = 0 (curves i) and a = 9 (curves 2) at U 1 /U 2 = 3/2; b — 
a = 9 at U 1 /U 2 = 6/5. 

corresponds to stability of flow for a = 0 due to the positive radicand in formula for c{k) 
for all k > 0. The flow with a ^ 0 loses stability for the sufficiently large values k. In the 
case considered /c* ~ 353. 

After analysis described above we conclude that long-wave perturbations (fc —)■ 0) 
are stable in a viscous fluid {a > 0). It gives grounds for the study of layered flows in 
long-wave approximation in details. 

4 Layered flows 

We consider the motion of an incompressible fluid in a Hele-Shaw cell. Its horizontal size 
L is much more than its depth H, i.e. H/L = e « 1. Further we assume 

t —)■ X e~^x, V ^ ev, fi ^ efi 

and omit all terms of order in equations ([3]). It corresponds to the long-wave approxi¬ 
mation [151 US] • As the result we have the model 

p{ut + (3{uu^ -F vUy)) +Px = -pu, Py = -gp, 

Ux + Vy = 0, Xf -|- uKx + VKy = 0; (13) 

p = p(x), /i = /i(x); = = 

Here |/ = 0 and y = H are the lower and upper boundaries in the ^/-direction. Note that 
stability of Hele-Shaw shear flows for weakly compressible barotropic fluid for /3 = 1 is 
studied in mi in the framework of hyperbolicity of integro-differential motion equations 
[T31 US). However for the long-wave models of incompressible inhomogeneous fluid the 
transition to the semi-Lagrangian variables does not lead to simplihcations of motion 
equations and does not apply here. 

We deal with the class of layered flows 


u = Ui{t,x), = const, {y e i = l,...,N), 











wherein 0 
form 


I/O < yi{t,x) < ... < yN{t,x) = 


pi (uit “h (3UiUix y “1“ 

V j=l / 

N 


H. For this class equations ffT3|l take the 

N 


9 hjx POx Pi^ii 

j=i+l 

N 


(14) 


X O 5 ^ ^ H.j ^ ^ Uihi , 

i=l i=l 

where hi = yi — yi-i is the depth of layer i with the density pi and the velocity Ui] Q is the 
constant flow rate through the vertical section of the cell; Pq is the pressure on the upper 
lid. For the derivation of the equations we use the kinematic condition on the interface. 
The similar models for unstratihed fluid {pi = p = const) are studied in [14] . 

System (na is reduced to evolutionary form using new dependent variables Sj = 
PiUi — PnUn- Then the motion equations are rewritten as follows 

(5( ' N N-l \ 

Sit + |( Piu] + Pnu^! + 2gpi ^ hj - 2gpN ^ hj + 2g hj \ = 

V i=i i=i i=i+i / X 


— —PiUi + PnUn, 


N N 

hit T i^^ihi^x O 5 ^ ^ hi hf, ^ ^ Uihi % 1 , . . . , N 1 . 

i=l i=l 

Depth and velocity of the last layer have the form 


N-l N-l N-l , 

hN = H-Y,h„ 

j=i j=i hj P] 

In some cases it is suitable to use a moving coordinate system. It has the average flow 
velocity U = Q/H. The corresponding change of variables is x' = x — Ut, u' = u — U, 
U = const. Then the hrst equation in (fT4|) takes the form 


^ 1 

+ (/3 — l)U^Uix + — hjx H- Pox = — — {ui + f/). 

Pi Pi Pi 

The primes are omitted. The other equations in (1141) do not vary. Further we consider 
modihcations of system (1T4|) and make numerical calculations of the interface position 
between two fluids with different physical properties. In particular it is discussed that in¬ 
terpretation of Saffman-Taylor instability can be given with the help of the one-dimension 
models. 


4.1 Two-layer flow of stratifled fluid 

Equations flT4|) for two-layer flow in the moving coordinate system have the form 
Piuit + Pi(/3mi + 'lU)uix + gpihix + gp2h2x + Pox = -Pi{ui + U), 

P2U2t + P2(^M2 + lU)u2x + POx = -P2{U2 + U), (15) 

hit + {uihi)x = 0, h2t + (^ 2 ^ 2 ) 3 ; = 0, 
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where 7 = /3 — 1. The depth and velocity of the second layer are given by formulas 


h2 = H — hi, U 2 =—{H — hi) ^uihi. 

We linearize the momentum equations in fllSp and obtain the modihed system 
PiMn + Pi7^“ix + 5'(Pi-P2)hix+Pox = -iJ.i{ui + U), 
p2U2t + P2IUU2X + POx = —^2(^2 + U). 

Using the change of variables mentioned above s = piUi — P 2 U 2 we obtain the motion 
equations in the evolutionary form 


St + (yUs + g{pi - p2)h)^ = {p2 - hi)U 


P2h + Pi{H — h) 
{H — h)pi + hp2 


(16) 


where 


ht {p^h)x 0 , 

h = hi, u = Ui = 


{H - h)s 


(H — h)pi + hp2 

It is not difficult to rewrite equations ffTBD in the form + AU 3 ; = F, where U = (s, h)"^ 
is the vector of unknown variables. The right part F and matrix A have the following 
form 

{p 2 h + pi{H - h))s 


F= - 


A = 


{H — h)pi + hp2 
lU g{pi - P 2 ) 


+ iP-2 — hl)U, oj , 

{H-h)hs 


V’ = 


{H — h) Pi + hp2 
Direct calculations show that the matrix A has the eigenvalues 


I f jj dip 


7 U 


dip\^ 


+ 4:g{pi- P 2 ) 


dip 

ds 


Correspondingly the characteristics of system flTHl) can be obtained from the solution of 
equation dx/dt = Ai, 2 - Due to 0 < h < H and ip^ = ip/s > 0 system flTHD is hyperbolic in 
the case of pi > p 2 - Standard methods developed for integrating the hyperbolic conserva¬ 
tion laws can be used for the numerical solution of equations flTB]) . We apply the simple 
and robust Nessyahu-Tadmor scheme for calculations of the interface position between 
two fluids with different physical properties. 

If we assume that the summands containing the velocity derivatives are equal to zero 
in equations flT^ then we have Darcy-like model 

g{pi - P2)hix + Pox = -pi{ui + U), pox = -P2{u2 + U), 

hit +{uihi)x, = t), U2 =-{H-hi)~^uihi. 

We can eliminate the pressure from these equations and specify the velocity ui as the 
function of layer height h = hi and its derivative hx- Now the system reduces to Burgers- 
like equation 

ht+{^{h,hx))^ = t), (17) 
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where 


$(h, hx) = uh = 


[H - h)h 


/i 2 - ^ii)U - g{pi - P 2 )K 


{H - h)pi + hp2 

If a fluid moves without gravity action [g = 0) or without stratihcation [pi = P 2 ) the 
equation obtained is replaced by the kinematic-wave model where $ does not depend on 
hx • 


Time 3 


Time 9 



We use parameters U = 1, if = 1 and /9 = 6/5 for calculations with respect to models 
(USD and f| 1 7p . We assume that initial interface position is as follows 


h 


t=Q 


1, X < Xq 

0, X > Xq. 


It is shown in Figs. [Hand [5] by the dashed line. We suppose additionally s(0,a;) = 0 for 
equation (ITSH . Results for model ([T7|) are presented in Fig. 0] at t = 3. It corresponds 
to displacing more viscous fluid (p 2 = 2) by less viscous one (pi = 1). The flow without 
gravitation [g = 0) is shown by curve 1. In this case we have the Saffman-Taylor instabil¬ 
ity, i.e. the interface position changes due to less viscosity of displacing fluid. Note that 
in the case of displacing more viscous fluid pi > p 2 the interface holds the initial position 
for all time a.t g = 0. Curve 2 corresponds to g = 980 and shows that positive stratihca¬ 
tion fades up the viscous huger velocity. Fig. [5] displays Rayleigh-Taylor instability (more 
viscous and denser huid displaces less viscous and less dense one due to gravitation) in the 
framework of model 0171) . Although displacing huid is more viscous the interface deviates 
from the initial position because the gravity prevails for the parameters given. 

Fig. |6] presents comparison of calculations of the interface position with respect to 
models flT6|) (curves 2) and f[T7|) (curves 1). Fig.[6l a corresponds to t = 9 and demonstrates 
change of the interface position without gravity action. It is not difficult to see that 
qualitative behavior of curves calculated with respect to diherent models coincides. If we 
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Time 9 



Figure 6 : The interface at ni = 1 , pi = 1 , 1 x 2 = 
2 corresponds to (fT6]) : a — ff = 0 , b — g = 980 


Time 3 jj 



change displacing and displaced fluids then instability does not develop and the interface 
does not move. 

There is the interface position for the gravity action in Fig. [ 6 l 6 at f = 3. The calcula¬ 
tion results coincide in the main flow region and differ in the vicinity of walls only. Note 
that flow structure for hyperbolic model 0161) has strong discontinuities adjacent to non¬ 
constant solution region. It is caused by the characteristics of equations ffT 6 |) are linearly 
degenerate at the initial data considered [19]. Truly the eigenvectors corresponding to the 
eigenvalues Xj (j = 1,2) of matrix A have the form = (Aj — ■ We construct 

the vectors = VAj. The gradient is calculated with respect to dependent variables s 
and h. Direct calculations show that scalar products ■ rj are equal to zero at s = 0 and 
h = 1 (h = 0). It explains the qualitative behavior of system ffTHj) solution in Fig. [ 6 l b. 

5 Conclusion 

The class of shear flow stability is studied for model ([3]). It describes Hele-Shaw flows of 
an inhomogeneous incompressible fluid taking into account dependence of viscosity and 
density on concentration. The latter holds along the trajectories. As a result of linearized 
equation solution the eigenvalues problem is formulated. In particular case this problem 
reduces to well-known Rayleigh equation. Analysis of two-layer flows of binary mixture 
with piecewise constant velocity, viscosity and density is provided (Fig. [T]). Boundary 
conditions at the interface ca and flTT]) are obtained. These formulas with equations 
(j9|) lead to dispersion relation ffT2|) . Unlike a number of other papers the factor = 6/5 
in front of convective terms is taken into account. This factor is generated at averaging 
through a cell gap. The dispersion curves are presented in Fig. [2] and demonstrate that 
viscosity stabilizes the long-wave perturbations. If the factor /? 7 ^ 1 then flow with slip 
line (contact discontinuity) without viscosity effect can be stable in a certain range of 
parameters. The stability region extends for the stratified flows (Fig. [3]). The growing 
perturbations exist for the wave numbers fc > /c*, where fc* is the critical number. Walls do 
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not influence on qualitative character of results in the framework of the scheme considered. 

With the help of equations of layered flows it is shown that the simplihed models are 
very suitable for the interface position modeling. The calculations with respect to these 
models demonstrate that the interface is not stable when more viscous fluid is displaced 
by less viscous one (at least for fluid without density stratihcation). The hyperbolicity of 
linearized model ffTHj) is proved for two-layer fluid at certain density relation. In the case 
of model based on Darcy law ffT7|) it is shown that it can be reduced to one Burgers-like 
equation. The kinematic-wave model (Hopf equation) can be obtained for homogeneous 
fluid with respect to density or without gravity action flow. The numerical calculations 
of the interface position using models ffTHD and ffTT)) display qualitative coincidence of the 
results. They differ in slower interface velocity close walls in the case of hyperbolic system 
ffT6|) only. 
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